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Abstract 

It has gone almost unquestioned that superexchange in the t — J (or Hubbard) model, and 
not phonons, is responsible for the unconventional ("d-wave") pairing symmetry of cuprate super- 
conductors. However a number of advanced numerical studies have not found superconductivity 
in the Hubbard (or t — J) model. On the other hand compelling experimental evidence for a 
strong electron-phonon interaction (EPI) has currently arrived. Here I briefly review some phonon- 
mediated unconventional pairing mechanisms. In particular the anisotropy of sound velocity makes 
the phonon-mediated attraction of electrons non-local in space providing unconventional Cooper 
pairs with a nonzero orbital momentum already in the framework of the conventional BCS theory 
with weak EPI. In the opposite limit of strong EPI rotational symmetry breaking appears as a 
result of a reduced Coulomb repulsion between unconventional bipolarons. Using the variational 
Monte-Carlo method we have found that a relatively weak finite-range EPI induces a d-wave BCS 
state also in doped Mott-Hubbard insulators or strongly-correlated metals. These results tell us 
that poorly screened EPI with conventional phonons is responsible for the unconventional pairing 
in cuprate superconductors. 
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I. INTRODUCTION 



It has been thought for a long while that Cooper pairs in the Bardeen-Cooper-Schrieffer 
(BCS) theory [l| with the conventional electron- (acoustic) phonon interaction are singlets 
and their wave function is isotropic (s-wave). This interaction has been thought to be local 
in space, so it could not lead to a higher angular momentum pairing. The pairing symmetry 
breaking is a many-body effect in accordance with a well-known quantum mechanics theo- 
rem, which states that the coordinate wave function of two particles does not become zero 
(or has no nodes) in the ground state Hence any superconductor should seem to be 
s-wave in the strong-coupling limit, where pairs are individual (e.g. bipolarons 3|, |J]) rather 
than overlapping Cooper pairs. 

Recently I have revised the symmetry of the superconducting state mediated by con- 
ventional acoustic phonons [5j. The sound speed anisotropy leads to a non-local attraction 
between carriers and unconventional Cooper pairs in the BCS layered superconductors in 
a wide range of carrier densities. The Bose-Einstein condensate (BEC) can also break the 
rotational symmetry due to a reduced Coulomb repulsion between unconventional small 
bipolarons bound by strong EPI. 

Earlier we have proposed that a strong departure of the cuprate superconductors from 
conventional BCS/ Fermi-liquids originates in the Frohlich EPI of the order of 1 eV [s], 1^, 7|, 
routinely neglected in the Hubbard U and t — J models of cuprate superconductors [8|. 
This interaction with c— axis polarized phonons is virtually unscreened because the upper 
limit for the out-of-plane plasmon frequency (< 200 cm^^ jo]) in cuprates is well below 
the characteristic frequency of optical phonons, c^o ~ 400 - 1000 cm Since screening is 
poor, the magnetic interaction remains small compared with the Frohlich EPI at any doping 
of cuprates (see also Ref . ) . Taking into account that the direct Coulomb repulsion is 
of the same order as the Frohlich EPI, we have proposed a so-called Coulomb-Frohlich 



model (CFM) of cuprate superconductors with the ground state in the 



small bipolarons or polaronic Coop er pairs (depending on doping) 0, [ill, 12 1, which can 
condense at high temperatures 13(]. More recently we have shown that even a weak long- 
range EPI combined with the Hubbard U provides sizable superconducting order in doped 
Mott-Hubbard insulators and/or strongly-correlated metals 14 1. 

Now compelling experimental evidence for a strong EPI has arrived from isotope effects 



brm of mobile 
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15l |. high resolution angle resolved photoemission spectroscopies (ARPES) [16|, a number 



of optical 



17| , neutron-scattering [l8|, |l9| and some other spectroscopies of cuprates |20| , in 



particular from recent pump-probe experiments [21|. These experimental observations and 
our theoretical findings tell us that EPl with conventional phonons is responsible for the 
unconventional pairing in cuprate superconductors. 



II. UNCONVENTIONAL COOPER PAIRS GLUED BY ACOUSTIC PHONONS 



It has gone unquestioned that the unconventional pairing requires unconventional 



electron-phonon interactions with specific optical phonons [22|, |23|, |24, l25l . I26l| . some- 
times combined with anti- ferromagnetic fluctuations [27] or vertex corrections 28|], or non- 
phononic mechanisms of pairing (e.g. superexchange pi), and a specific shape of the Fermi 
surface. Here I show that even conventional acoustic phonons can bound carriers into un- 
conventional Cooper pairs due to the sound-speed anisotropy in layered crystals. 

In the framework of the BCS theory the symmetry of the order parameter A(k) and the 
critical temperature, T^, are found by solving the linearised "master" equation, 

k 



A(k) = -5^l^(k,k')^tanh 



2kBT, 



(1) 



The interaction l^(k, k') comprises the attraction, —Vph{q), mediated by phonons, and the 
Coulomb repulsion, K(q) as V(k,k') = -Vph{q)Q{u;D - \^k\)0{u;D - l^k'l) + V,{q)Q{Up - 
|^k|)0(^^p — l^k'l), where Vph{q) = C^/NMcf is the square of the matrix element of the 
deformation potential, divided by the square of the acoustic phonon frequency, cjq = Qg, 
Q is the sound speed, M is the ion mass, is the number of unit cells in the crystal, and 

is the electron energy relative to the Fermi energy. The magnitude of C is roughly the 
electron bandwidth in rigid metallic or semiconducting lattices. The electron momentum 
transfer q = k — k' or its in-plane component has the magnitude q = 2^^'^kp[l — cosipY^"^ for 
the spherical or cylindrical Fermi surface, respectively, where ip is the angle between k and 
k', and hkp is the Fermi momentum. Theta functions account for a difference in frequency 
scales of the electron-phonon interaction, ujd, and the Coulomb repulsion, Up ^ ud, where 
uj£, and ujp are the maximum phonon and plasmon energies, respectively. 

If one neglects anisotropic effects, replacing Vph{q) and Vc(q) by their Fermi-surface 
averages, Vph{q) =^ Vph, K;(q) =^ Vc, then there is only an s-wave solution of Eq.([l]), 



As, independent of k. The sound speed anisotropy actually changes the symmetry of the 
BCS state. While q is a constant in the isotropic medium, it depends on the direction 
of q in any crystal. The anisotropy is particulary large in layered crystals like cuprate 
superconductors. As an example, the measured velocity of longitudinal ultrasonic waves 
along a — b plane, C||=4370 ms^^ is almost twice larger than that along c axis, c_l=2670 ms~^ 



in Bi2Sr2CaCu20 



8+y 



291]. It makes Vph{q) anisotropic. 



q2 

where a = (cy — cjj/c^ is the anisotropy coefficient, which is about 2 in cuprates. The 
corresponding real-space potential is non-local, 

5(r) a 



(3) 



d 47r(l + a)i/V3j ' 

falling as at the distance r ^ d between two carriers in the plane, where Vph = C"^ /Mc\. 
Also the Coulomb repulsion is q dependent, K(q) = 47re^/\^eo(5'^ + qT)- framework 
of the random phase approximation the inverse screening radius squared is found as = 
87re^A^(0)/Veo with the density of states (per spin), iV(0), at the Fermi surface. Here d is 
the inter-layer distance and eg is the {in-plane) static dielectric constant of the host cuprate 
lattice of the volume V . 

Solving the master equation ([T]) with 2D electron spectrum one can expand A(k) = 
Y,m ^rn exp(zm0) and VphA^) = Y.m ^p/i,c(g±, m) exp[zm(0 - 0')] in series of the eigenfunc- 
tions of the c-axis component of the orbital angular momentum, where and 0' are polar 
angles of the in-plane momenta, ky and kj|, respectively. The critical temperature of an 
m-pairing channel (m = 0, ±1, ±2, ...) is found as 

T,m = l.UcoD exp [-1/(A™ - fj.*J] , (4) 

where /xj!^ = + fJ-m^^i^p/^o)]- Here and are the phonon-mediated attraction 

and the Coulomb pseudopotential in the m-pairing channel, given respectively by 



A,. _ ^ a n dx[x + 1 - ^x{x + 2)]- 

A 2^ Jq Vx + 2 



and 

/im _ V7 r ^^[^ + /5 + 1 - ^/{x + l3){x + (3 + 2)\ 



l^c 2 Jo ^x{x + f3){x + f3 + 2) 



(6) 
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where A = N(S))C'^ /NMc\, 7 = 71"^ /2<Pkl{l + a), 7 = 7(1 + a), = Ae^<fN{<d)/'KVeo, 
and [3 = qg/2kp (note that A, /ic, and do not depend on the carrier density since A^(0) is 
roughly constant in the quasi- two dimensional Fermi gas). 

The effective attraction of two electrons in the Cooper pair with non-zero orbital mo- 
mentum turns out finite at any finite anisotropy, a ^ 0, but numerically smaller than in the 
s-channel, as seen from its analytical expression for s-wave, m = 0, p-wave, m = 1, d-wave, 
m = 2, and higher orbital momentum pair states, obtained by integrating in Eq.Q. The 
Coulomb repulsion turns out much smaller in the unconventional pairing states than in the 
conventional s-wave state, which is seen from the analytical expression for fj,m, Eq.(l6]) and 
from Fig.l. 

Using the simplest parabolic approximation for the 2D-electron energy spectrum we can 
draw some conclusions on the carrier-density evolution of the order-parameter symmetry. 
Within this approximation, kp = 27Tdn and A^(0) = m*V/27rdh'^, where n = 2x/Q is the 
carrier density, m* is the effective mass, and x is the doping level as in La2--xSra.Cu04 
with the unit cell volume Q. The ratio of the parameters (3 = m*e'^Q/27rh'^d'^eox and 7 = 
TrQ/Sd'^x ~ 0.044/x is independent of the carrier density, f3/'y = 4m*e^c//7r^/i^eo, which is 
approximately 5 for the values of = Arrie and eg = 10. Fixing the value of the EPI 
constant at A = fic/^'^ (which corresponds to the weak-coupling BCS regime with A ~ 0.1 
since fic is of the order of 1) and taking ln(co'p/ti;£)) = 3, we draw the anisotropy- doping 
phase diagram. Fig. 2, with the critical lines for s, p and d order parameters, defined by 
Am = l^m- The state with the lowest magnitude of the anisotropy, a/(l + a)^/^, is physically 
realized since it has the highest Tc. At substantial doping the screening length becomes larger 
than the typical wavelength of electrons, (3 ^ 0, so that the s-wave state is the ground state 
at a large number of carriers per unit cell for any anisotropy. On the contrary, the Coulomb 
repulsion is reduced to the local interaction at a low doping, /3 ^ 00, and rf-wave Cooper 
pairs are the ground state even at very low value of the anisotropy, Fig.2. Interestingly, s— 
and (i-states turn out degenerate at some intermediate value of doping, x = Xc- Hence there 
is a quantum phase transition with increasing doping from d— to s-superconducting state, if 
a > Oc, and from d— to the normal state and then to the s-wave superconductor, if a < Oc, 
Fig.2. 
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III. BREAKDOWN OF ROTATIONAL SYMMETRY IN THE STRONG- 
COUPLING LIMIT 



In the strong-coupling regime, A > 1, the pairing is individual [3|, in contrast with the 
collective Cooper pairin g. W hile BEC of individual bipolarons can break the symmetry 



on a discreet lattice 



30 



3l| . I have proposed a symmetry breaking mechanism [5], which 



works even in a continuum model, where the ground state, it would seem, be s-wave to 
satisfy the theorem. The unscreened Frohlich EPI in layered ionic lattices like cuprates has 
been suggested by us as the key for pairing [s!]. Acting alone it cannot overcome the direct 
Coulomb repulsion, but almost nullifies it since eg ^ 1. That allows the weaker deformation 
potential, Eq.([2]), to bind carriers into real-space bipolarons, if A ^ 0.5 While its local 
part is negated by the strong on-site repulsion U, the non-local tail provides bound pairs 
of different symmetries with the binding energies > Ap > A^ > ... in agreement with 
the theorem. However, there is the residual Coulomb repulsion between bipolarons, Vc{R)- 
Since bipolarons have a finite extension, ^, there are corrections to the Coulomb law. The 
bipolaron has no dipole moment, hence the most important correction at large distances 
between two bipolarons, ^, comes from the charge-quadrupole interaction, 

. 4e'i±^. (7) 

where ?7 is a number of the order of 1, and plus/minus signs correspond to bipolarons in the 
same or different planes, respectively. The dielectric screening, eo is highly anisotropic in 
cuprates, where the in-plane dielectric constant, eoii, is much larger then the out-of-plane one. 



321 ]. Hence the inter-plane repulsion provides the major contribution to the condensation 
energy. Since oc 1/A, the repulsion of unconventional bipolarons with smaller binding 
energies, A^, Ap < A^, is reduced compared with the repulsion of s-wave bipolarons. As 
a result, with increasing carrier density we anticipate a transition from BEC of s-wave 
bipolarons to BEC of more extended p— and d-wave real-space pairs in the strong-coupling 
limit. 



IV. CONCLUSIONS 



Several authors 8| have remarked that superexchange, and not phonons is responsible 
for the symmetry breaking in unconventional superconductors like doped cuprates. Here I 



arrive at the opposite conclusion. Indeed, superexchange interaction, J, is proportional to 
the electron hopping integral, t, squared divided by the on-site Coulomb repulsion (Hubbard 
U), J = 4t'^/U, estimated as J ~ 0.15 eV in cuprates jsj. This should be compared with the 
acoustic-phonon pairing interaction, Vph, which is roughly the Fermi energy, Vph ~ Ep ^ 4t 
in a metal, or the bandwidth squared divided by the ion — ion interaction energy of the order 
of the nearest-neighbour Coulomb repulsion, Mcf « 14 in a doped insulator. The small ratio 
of two interactions, J /Vph ~ t/U ^ 1, or J /Vph ~ Vc/U <^ 1 and the giant sound-speed 



anisotropy 



291] favor conventional EPI as the origin of the unconventional pairing both in 



underdoped cuprates, where the pairing is individual [3j, and also in overdoped samples 
apparently with polaronic Cooper pairs, which can coexist with bipolarons |l2l ]. 

Moreover recent studies by Aimi and Imada 3^ of the Hubbard model, using a sign- 
problem-free variational Monte Carlo (VMC) algorithm, have shown that previous approx- 
imations overestimated the normal state energy and therefore overestimated the condensa- 
tion energy by several orders of magnitude, so that the Hubbard model does not account for 
high-temperature superconductivity. This remarkable result is in line with earlier numerical 
studies using the auxiliary- field quantum (AFQMC) 3^ and constrained-path (CPMC) 35] 
Monte-Carlo methods, none of which found superconductivity in the Hubbard model. On 
the other hand using VMC method we have found that even a relatively weak finite-range 
EPI induces the d-wave superconducting state of doped Mott-Hubbard insulators and/or 
strongly-correlated metals with a sufficient condensation energy jl^ . 

I conclude that the finite-range electron-phonon interaction is the key to the high and a 
higher temperature superconductivity. 

I would like to thank Tom Hardy, Jim Hague, Victor Kabanov, Pavel Kornilovitch and 
John Samson for long-standing collaboration and illuminating discussions. The work was 
supported by EPSRC (UK) (grant numbers EP/C518365/1 and EP/D07777X/1). 
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FIG. 1: The Coulomb repulsion, firm as a function of the ratio of the electron wavelength to the 
screening length squared (/? = q'g/2kp), and the electron-phonon coupling constant, A^, as a func- 
tion of the ratio of the electron wavelength to the inter-plane distance squared, 7 = n'^ /2cPkp{l+a) 
for a = 4 (inset) in s,p and d pairing channels. Here fi'^ = ficj- (Reprinted with permission from 
Ref.^]. ©2008 by the American Physical Society.) 
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FIG. 2: Critical sound-speed anisotropy, a/(l + a)^/^ = (c| — c^)/c|| C_|_ , 3fS 3j function of doping, 
for A = /ic/12 (solid lines correspond to d and s states, and dashed line to p-state). With 
increasing carrier density there is a quantum phase transition at x = Xc from a d-wave to an s- 
wave superconductor, when a > Oc, and two quantum phase transitions from d-wave to the normal 
state and from the normal state to the s-wave state when a < Uc- (Reprinted with permission 
from Ref-Q]. ©2008 by the American Physical Society.) 
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